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Abstract. Using wavelet analysis approach, we can de- 
rive a measure of the disorder content of solar activity, 
following the temporal evolution of the so-called wavelet 
entropy. The interesting feature of this parameter is its 
ability to extract a dynamical complexity information, 
in terms of frequency distribution of the energy content, 
avoiding restrictions, common in the nonlinear dynamics 
theory, such as stationarity. The analysis is performed 
on the monthly time series of sunspot numbers. From 
the time behaviour of the wavelet entropy we found a 
clear increase in the disorder content of solar activity for 
the current 23 th solar cycle. This result suggests general 
low accuracies for current solar cycle prediction methods. 
Moreover, we pointed out a possible connection between 
wavelet entropy behaviour and solar excursion phases of 
solar dipole. 
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1. Introduction 

Multiscale analysis, based on the wavelet approach, has 
been successfully used for many different physical appli- 
cations, including glo bal so lar activity (Ochadlick et al. 
Lawrence et al. 1995). Intermittence features have 



1992 



been well analysed through the use of a continuous wavelet 
transforms on the monthly sunspot numbers time series. 
More recently, different authors used longer daily sunspot 
numbers time series to investigate more accurately new 



subtle periodicities and their evolutions (Frick et al. 1997; 
Ballester et al. |1999| ; Fligge et al. |1999| Sello |2000[ ). A re- 
cent work on wavelet applications to biomedical signals, 
shows the utility of a new parameter, derived from wavelet 
formalism, to quantify the degree of disorder in the time 
series. This quantity, called wavelet entropy, essentially 
gives an information on the scales extension involved in 
the energy distribution, i.e. a dynamical complexity pa- 



rameter (Quian Quiroga et al. 199£). The main aim of the 



present work is the application of the wavelet entropy ap- 
proach to the solar activity, through the wavelet analysis 
of the monthly mean sunspot numbers. 



2. Wavelet entropy 

Fourier analysis is an adequate tool for detecting and 
quantifying constant periodic fluctuations in time series. 
For intermittent and transient multiscale phenomena, the 
wavelet transform is able to detect time evolutions of the 
frequency distribution. The continuous wavelet transform 
represents an optimal localized decomposition of time se- 
ries, x(t), as a function of both time t and frequency (scale) 
a, from a convolution integral: 
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where ip is called an analysing wavelet if it verifies the 
following admissibility condition: 



where 
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is the related Fourier transform. In the definition, a and 
r denote the dilation (scale factor) and translation (time 
shift parameter) , respectively. We define the local wavelet 
spectrum: 
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where fco denotes the peak frequency of the analysing 
wavelet ip. From the local wavelet spectrum we can de- 
rive a mean or global wavelet spectrum, P UJ (k): 



dt P u (k,t) 



(5) 



The relationship between the ordinary Fourier spectrum 
Pf{oj) and the mean wavelet spectrum P u (k) is given by: 
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indicating that the mean wavelet spectrum is the aver- 
age of the Fourier spectrum weighted by the square of the 



2 



S. Sello: Wavelet entropy as a measure of solar cycle complexity 



Fourier transform of the analysing wavelet ip shifted at 
frequency k. Here we used the family of complex analyz- 
ing wavelets consisting of a plane wave modulated by a 
Gaussian, called Morlet wavelet, (Torrence et al. 1998) ) : 



n -l/A e i U0Ve -V>/2 
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where loq is the non dimensional frequency here taken to 
be equal to 6 in order to satisfy the admissibility condition, 
Eq.(2). Following Quian Quiroga ct al. (1999), we define 
a wavelet entropy as a function of time: 



WS t = - ^2,Pt,klo9i{pt,k) 



where: 
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is the energy probability distribution for each scale level. 
From the definition, follows that an ordered activity corre- 
sponds to a narrow frequency distribution of energy, with 
low wavelet entropy, and a random activity corresponds to 
a broad frequency distribution, with high wavelet entropy. 
Of course, higher values for wavelet entropy means higher 
dynamical complexity, higher irregular behaviour, lower 
predictability. The application of the wavelet entropy is 
optimal for non-stationary signals. 

3. Solar activity: sunspot numbers 

We considered here the record of solar activity 
given by the monthly me an number of sunspots 
(URL: [http: //www.oma.bc/KSB 



from SIDC archive 



ORB/SIDC/indcx.html ), covering the time interval: 1749- 
2000.04 and consisting of 3013 observations. Using a 
proper wavelet analysis we found the results displayed in 
Fig.l. Following the notation described in Sello ( 2000| ), 
the upper part shows the original time series in its natu- 
ral units. Time is here expressed in years. The central part 
shows the amplitudes of the wavelet local power spectrum 
in terms of an arbitrary greyscale contour map. White 
higher values are strong energetic contributions to power 
spectrum, while black lower values are weak energetic con- 
tributions. Horizontal time axis corresponds to the axis of 
time series and vertical scale (frequency) axis is, for conve- 
nience, expressed in log values of cycles per year -1 . Thus 
the range analyzed is between 148 years (value -5) and 
134 days (value 1). The right part shows the mean global 
wavelet power spectrum (solid line, SM) obtained with a 
time integration of the local wavelet power spectrum, and 
the 5% significance level using a red noise autoregressive 
lag-1 background spectrum with a = 0.93 (dotted line, 
SC) More precisely, we assumed here that different re- 
alizations of the considered stochastic process, given by 
sunspot activity, will be randomly distributed about a red 



Fig. 1. Wavelet analysis of monthly mean sunspot num- 
bers - right: mean global wavelet power spectrum (solid 
line, SM); and 5% significance level (dotted line, SC) 



noise background spectrum (i.e. increasing power with de- 
creasing frequency) , which is used to verify a null hypoth- 
esis for the significance of a given peak in the wavelet 
power spectrum. In particular, if a peak in the wavelet 
power spectrum is significantly above the assumed red 
noise background spectrum, then it can be considered a 
real feature of the original signal with a given confidence 
level. This statement appears well justified from different 
applications to many complex physical signals, e.g. geo- 
physical data. In order to obtain a simple model for red 
noise we considered a univariate lag-1 Markov process: 
x n — otXn-i + z n , where a is the assumed lag-1 autocor- 
relation, and z n is a Gaussian process. Note that: a = 
correspond s to a white noise, (for more details see: Tor- 
rence et al. 1998 ). From the figure, it is possible to distin- 
guish the main region corresponding to the well defined 
Schwabe (~ 11 years) cycle, with a dominant energetic 
contribution above the red noise background spectrum. 
In particular, the complex evolution of the cycle, both for 
scales and amplitudes, is well evident. The irregular finite 
extension in the map represents a continuous interaction of 
frequencies which generate complex dynamical behaviour. 
In order to better quantify the evolution of complexity, or 
disorder for the solar cycle time series, we computed the 
time evolution of the wavelet entropy. An interesting fea- 
ture of this approach is that the usual assumed hypothesis 
of stationarity of time series is here not necessary. 

Fig. 2 (lower part) shows the time evolution of the 
wavelet entropy, Eq.(8). For comparison, an ordered de- 
terministic (harmonic) signal gives the numerical result 
displayed by the dashed line. It is interesting to note that 
the degree of disorder reaches maximum local values, often 
in phase with solar maxima, where the magnetic activity 
results higher, but with many important exceptions (see, 
for example, cycles 20 and 22). The recorded maximum 
entropy value has been reached during the maximum of 
solar cycle 4 (1789). Moreover, the entropy results quite 
low for long periods, including solar cycles from 5 to 7 
(1798-1828) and, in particular, for solar cycles from 9 to 
14 (1848-1912). For the last three cycles, 21, 22 and 23, we 
found a clear average increase of the wavelet entropy, sug- 
gesting a more complex dynamics, with a higher level of 
disorder due to a broad frequency energy distribution, con- 
nected to evolution of the global magnetic activity. This 
property forces a reduction in the accuracy level reachable 
by prediction methods, which are based on some strong 
deterministic hypothesis. For nonlinear chaotic dynam- 
ics methods this means, at least, a low er value for Lya- 
punov's predictability time, (Sello 1999 ). Furthermore, we 
note that the level of wavelet entropy can help to quantify 
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Fig. 2. Monthly sunspot numbers (up) and wavelet en- 
tropy (down) time evolutions. The wavelet entropy for an 
ordered harmonic signal, is also displayed (dashed line). 

Fig. 3. Wavelet analysis of the wavelet entropy - right: 
global power spectrum (solid line, POW); and 5% signifi- 
cance level (dotted line, SL). 



all increase of the wavelet entropy for the last cycles. The 
suggested correlation between wavelet entropy and solar 
excursion phases is here limited to a pure qualitative view- 
point. Further investigations, with a more quantitavive 
evaluation based on geomagnetic indices, are thus needed 
in order to support a reliable strict link between the de- 
gree of disorder in the sunspot activity, as measured by 
wavelet entropy, and the solar dipole mechanism. 



the intermittence degree in solar activity caused by all the 
stochastic sources and, in particular, by stochastic fluctu- 
ations in meridional circulation of plasma, which are able 
to influence significant unpredictable variation in ampli- 
tude and phase of the solar cycle , as recently evidenced 
by Dikpati & Charbonneau ( |1999|) . Finally, we point out 
that the new information, added by the wavelet entropy, 
it is not useful for sunspot predictions. 

In order to investigate more accurately the periodici- 
ties which drive the evolution of the wavelet entropy, we 
performed a wavelet analysis on the above derived time 
series. For this analysis we used a wavelet software pro- 
vided by C. Torrence and G. Compo, and av ailable at 



URL: http: / /paos. colorado.edu/research/waveleti 



Fig. 3 depicts the wavelet map with 5% confidence level 
regions. Here the scale axis shows the periods in log val- 
ues. The u-shaped line corresponds to the so-called Cone 
of Influence, i.e. the region of the wavelet spectrum in 
which edge effects, due to the fi nite-le ngth of signal, be- 
come important. (Torrence et al. 1998| ). The global power 
spectrum (right part, POW) indicates that there is mainly 
one significant persistent periodicity related to the wavelet 
entropy, and localized around 40 years, a possible period 
doubling value of the main Schwabe cycle (Lawrence et al. 
1995). Other detected lower periodicities are more local- 



ized and energetically weaker, below the 5% significance 
level of the assumed red-noise background spectrum (right 
part, SL). The time evolution of the wavelet entropy sug- 
gests also the existence of other longer modulations, not 
resolvable in the available time series. 

It is interesting to point out possible relations exist- 
ing between wavelet entropy behaviour and solar excur- 
sion phases, recently investigated by Mursula et al. ( 1998| ). 
The excursion phases of the solar dipole are linked to the 
strong half solar rotation periodicity and arise when the 
heliospheric current sheet is flat and tilted. In particular, 
the authors found a clear connection between sunspot ac- 
tivity and the solar dipole tilt. More precisely, the overall 
occurrence of excursions remained very weak during solar 
cycles from 9 to 14, where the solar dipole moment was 
lower with a high stability of the coronal evolution, and 
this is coherent with the lowest values of wavelet entropy. 
Instead, during the most recent cycles (15-22) we detected 
largest excursions of the solar dipole, with a different cy- 
cle distribution, larger instabilities and a more disturbed 
heliosheet. This behaviour is again coherent with the over- 



4. Conclusions 

From wavelet analysis approach, we considered a mea- 
sure of the disorder content of solar activity, following 
the temporal evolution of the so-called wavelet entropy. 
The analysis is performed on the monthly time series of 
sunspot numbers and shows an interesting multiscale evo- 
lution of the wavelet entropy. In particular, we found a 
clear increase in the disorder content of solar activity for 
the current 23" 1 solar cycle. This result suggests general 
low accuracies reachable for current solar cycle prediction 
methods, which are based on some strong deterministic 
hypothesis. Finally, we suggested a possible strict corre- 
lation between the wavelet entropy qualitative behaviour 
and the solar excursion phases of solar dipole, linked to 
coronal complexity evolution. 
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